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In this paper, we prove that the ground state of the Kondo model with large spin is nondegenerate, 
apart from a SU(2) spin degeneracy in the case of half filling. The ground state spin is found for 
the system, and the energy level orderings are discussed. Finally, the existence of ferrimagnetism in 
some cases is proved. 
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The Kondo models, or single and lattice impurity mod- 
els, are one of the_most challenging subjects in strongly 
correlated systemsErO. Recent developments of reflection 
positivity technique in the spin space make it possible to 
establish some rigorous results for the half-filled strongly 
correlated electron systemsQ □. Theory of reflection pos- 
itivity in the spin space for the single- and multi-channel 
Kondo models with spin 1/2 was developed recently and 
a series of rigoroua-jresults on the ground state prop- 
erties were proved™!. However some materials are de- 
scribed by the Kondo models with large spin, such as 
{Lai^ x X x )MnOz with X = Ba, Ca, Sr etc. and the lo- 
calised spin s — 3/2E3. Due to more degrees of freedom 
in the case of large spin than in the case of spin 1/2, usu- 
ally it is very hard to extract rigorous results for those 
systems. As a generalisation of the theory for the Kondo 
model with spin 1/2, we will investigate the Kondo model 
with large spin and provide some rigorous results on the 
ground state of the Kondo model in the case of half fill- 
ing. 

Let us first write down the Hamiltonian we will inves- 
tigate: 

2 l 

H = E ^3 c Xa c 3,<y + Ui(m,i - -)(n hl - -) 

(ij>SA,CT ieA 

+ S JiSi ■ S c i + ^ KijSi ■ Sj, 

i£A ii£A d 



(1) 



where c[ a and Cj )<7 are the creation and annihilation op- 
erators for the conduction electron (c-electron) at site 
% with spin er(=f,J,) and nj )(T = c\ a Ci tC! . Si is the 
localised spin operator with spin s, at site i. S C i = 

J2a,a> c l,o-°V,<>-'ei,a-'/ 2 ! and a a (a = x, y, z) are the Pauli 
matrices. The model is defined on a bipartite lattice A 
with the site numbers Na and Nb of the two sublattices 
A and B. tij = tji is possibly non-zero only when i and j 
belong to two different sublattices. The lattice A is con- 
nected by the hopping terms {tij}, which implies that for 
any two sites k and / on A, we can always find a sequence 
{(Mi), (h,i2), • • ■ , (i n , I)} such that t kil t ili2 ■ --t ln i / 0. 
Ad is the distribution of localised spins on A. If we re- 
gard the same site for conduction electron and localised 
spin as two independent sites in a generalised lattice, the 
site i for conduction electron and the site i for localised 



spin belong to the same sublattice when J; < 0, and to 
two different sublattice when J ; > 0. In the case the 
generalised lattice can be still regarded as a generalised 
bipartite lattice. Assume KijJiJj < if i and j belong 
to the same sublattice, and K^^Jj > if i and j belong 
to two different sublattices. When Ji and Jj are constant 
or have the same sign, Kij < if i and j belong to the 
same sublattice and > if i and j belong to two dif- 
ferent sublattices as in the antiferromagnetic Heisenberg 
model on a bipartite lattice. This condition guarantees 
that the last term in Eq.(l) does not introduce any frus- 
tration for the generalised bipartite lattice. Physically, 
it is completely absent of frustration in this case. The 
model is reduced to the Hubbard model when A^ = 0, 
i.e., there is no magnetic impurity or localised spin in the 
systems. The structure function of the bipartite lattice 
A is defined as e(i) = 1 if i S A and — 1 if i G B. 

The main results are summarised as follows: 
Theorem: Assume the model in Eq.(l) with and 
is defined on a connected bipartite lattice A with 
sublattice sites Na and Nb and with the distribution 
of localised spins Ad- All Ui > and Ji and the 
number of conduction electrons is iV e = N A = Na + Nb ■ 
The signs of Kij and Ji satisfy the previously stated con- 
dition. Denote the lowest energy state in the subspace 
decomposed by the z-component of the total spin Sf ot by 
|*). 

i) . The state |^) is non-degenerate in each subspace of 
Sf 0t . The ground state of Eq. (1) is unique apart from a 
(2S° ot + l)-fold spin degeneracy. 

ii) . The total spin 5 tot in the lowest energy state is 



St, 



c0 

\s:. 



if \S? 



< c0 



otherwise 



where 



Si-^-re(i) 



Ji 



iSA iSAd 

iii). The spin-spin correlation functions obey, 
(tf|S+ • S".|*> = e(i)e(j)C ir , 
(*\Sf-Sj\*)=e(i)e(j)F ij ; 

(*|S+.S7|*) = -c(i)e0")|7rG ! y, 



(2) 



(3) 



(4) 
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where Cy, Fij and Gij > if Ui > 0, and > if all 
C/i > 0. 

iv). When all Ui = 0, (at least one of) the ground state(s) 
(if degenerate) has the total spin as that in Eq. (0). 

Before we present the proof, several remarks or corol- 
laries are made: 

1). In the case of all Ji > or < and s, = s, suppose 
that NacI spins are on the sublattice A, and Nsd spins 
on the sublattice B, the total spin is 



c0 _ 
D tot — 



(N A - N B ) - -irs(N Ad - N Bd ) 



J 



(5) 



In the Kondo lattice case, Na 
and the total spin is 



N Ad and N E 



N 



c-0 
°tot 



1 J 

2 L7j 



\n a 



Bd; 



(6) 



When all Sjr-j= 1/2, we recover the result for the case 
of spin l/2Ertl The ground state is a singlet only when 
N A = N B , or s = 1/2 and J > 0. When N A # and 
(A^t — Nb)/{Na + Nb) ^ when the system becomes 
sufficiently large, we obtain a state with ferromagnetic 
long-range order. 

2). Theorem (hi) indicates that strong antiferromag- 
netic correlations exist between the conduction electrons 
or localised spins. A direct corollary is that the antifcr- 
romagnetic correlation is always stronger than the ferro- 
magnetic correlation. For example on a cubic lattice, 



(S+.S Q )=max{(S+.S: 9 )}, 



(7) 



where S + = 



and Q — (tt, 7r, • • •). 



If the state possesses ferromagnetic long-range order, it 
must also possess antiferromagnetic long-range order. In 
other words, the ferromagnetism in the case of Na ^ Nb 
is, strictly speaking, the ferrimagnetism. 

3). The energy level orderings can be obtained from 
Theorems (i-ii) and the SU(2) symmetry with a varia- 
tional principle: denote E(S% ot ) the lowest energy state 
with Sf ot . As \^(S z ot )) is the lowest energy state 
with S' ot , we can construct an eigenstate ($Z igA S~ + 
£ieA d K)\^(Stot)) (suppose S z tot > 1) due to the spin 
SU(2) symmetry. This state has its z-component of the 
total spin S z ot — 1 and the eigenvalue E(S z ot ). Meanwhile 
it has the same total spin as in |^(5j ot )). In the varia- 
tional principle the lowest energy state in the subspace 
of S z ot — 1 should not be higher than E(S^ ot ). If the low- 
est energy in the subspace S z ot — 1 is equal to the lowest 
energy in the subspace S* ot , the lowest energy state or 
one of the states if degenerate with S^ ot — 1 must has the 
same total spin in the state with S^ ot . Since the lowest 
energy states in each subspace are non-degenerate when 
all Ui > and we assume S* ot > 0, we have 



E ( S tot) = E {- S tot)\ 

E{S z tot ) < E(S? ot + 1) if SI 
E{S* tot ) = E(S° ot ) if 0< SI 



£- °tot ! 
^ Stot ■ 



(8) 
(9) 
(10) 



For S* ot < 0, the energy level ordering can be obtained 
from Eqs. £8rl0). This is similar to that of the Hciscn- 
berg modelt3. It is worth of mentioning that in the case 
of Ui > we have to use "lower or equal" instead of 
"lower than" if we could not determine whether the low- 
est energy state in each subspace is non-degenerate or 
not 

4) . In the one-dimensional chain of the Kondo lattice 
with all Ji > or < 0, Theorems (ii) and (hi) on the to- 
tal spin and spin-spin correlation functions are the same 
as those in a two-chain spin-ladder system. This coin- 
cides with the analysis on thc_L'csemblance of these two 
systems by White and AfflccktLil 

5) . When the system has no localised spin or magnetic 
impurity, NAd — Nsd — and the model in Eq.(l) is 
reduced to a Hubbard model. Theorem (i-iv) hold for 
the Hubbard model. In the case, one can regard the spin 
formula in Eq. (2) as a. generalisation of Lieb's theorem 
for the Hubbard modelcl. The conditions of even number 
of the lattice sites and 5f ot = in Lieb's theorem are re- 
moved. As a corollary, the energy level ordering for the 
Hubbard model in the case of half filling is the same as 
shown in Eqs. (8-10). 

The main steps to prove this theorem are 1). to ex- 
press the spin operator in a multi-fermion representation; 
2). to define the generalised bipartite lattice after we 
consider the localised spins; 3) to introduce a complete 
and orthonormal set of basis for the system; 4) to prove 
the positive definiteness oiiscmidcfinitencss of the ground 
state on the chosen basisllj; and 5) to prove (hi) and (iv) 
by utilising the positive definiteness of the ground state 
and the known theorems on a positive definite state. 

Except for the reflection positivity_approach in the 
spin space as in the case of spin 1/2Q'lI, the key trick 
used in this paper is to express the localised spin op- 
erator Si as the summation of 2si spin 1/2 operators 
in a multi-fermion representation (d-fermion) and these 
2si spins are coupled ferromagnetically. Completely fer- 
romagnetic coupling could be realized by introducing a 
strong ferromagnetic coupling limit between these spins 
in the Hamiltonian. The spin operator is written as 



2s 4 

E 



(ii) 



with the restriction that each site ion is singly occupied 
by d-fermions. The effective Hamiltonian in terms of d- 
fermions is 

(ij)e.A ieA 
- E A »[(E A^^di^f - s t (s t + 1)], (12) 
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where all A, are positive^. We shall take all Xi — > oo 
at the last step (as a matter of fact our theorem is true 
for any finite Xi), and this guarantees that any deviation 
from 



E 4 a J 



= Si( Si + 1) (13) 



will lead to the divergence of the eigenvalues of energy in 
Eq. (12). Eq. (13) is an alternative expression of 

S?|*) = *( ai + l)|*) (14) 

in the multi-fermion representation. In the limit of 
Xi — > +oo, H' in Eq. (12) is equivalent to H in Eq. 
(I)- 

The definition of the generalised bipartite lattice is 
very crucial for utilising the reflection positivity tech- 
nique in the spin space. The generalised bipartite lattice 
for H' is defined as follows: 1). The sublattice A (B) 
for c-electrons belongs to the generalised sublattice A 
(B); 2). The site icti for d-fermions belongs to the gen- 
eralised sublattice A (B) if i belongs to the sublattice A 
(B) when J, < and to the generalised sublattice B (A) 
when Ji > 0. At the same site i all icti with different a* 
belong to the same generalised sublattice. Additionally, 
the sites of the generalised lattice which are connected 
by Kij belongs to the same sublattice when < and 
to two different sublattices when Kij > according to 
the condition on the signs of Ji and K^ . 

According to this generalised bipartite lattice, the 
transformation operator T forr-the partial particle-hole 
transformation is introduced astil 



T=n(*T-*(i)4r) n n (<* 



i<EA 



Hat 1 



(15) 



Under this transformation, we have 



Tc lT Tt 


= (- 


-r 




(16) 


Tci+Tt 


= (- 


-i; 




(17) 




= (- 


-i; 




(18) 




= (- 


-r 




(19) 



where N t is the total number of the generalised bipartite 
lattice sites. 

The set of basis for W in the case of half filling 
we choose is {T|^) ® \<f>p}}- {IO} is a rea lj com- 
plete and orthonormal set of basis for No c-electrons 
and d-fermions with spin a, which is expressed as (\f a — 

Ihea C L U. ieai 4a ie a|0) ( the OTder of c and d operators 
will cause an undetermined phase — 1 , which does not af- 
fect our final result). The single occupancy of d-fermion 
at site ion is realized by choosing all {a, /3} such that 



n-ia^l) ® |<^} = n io!il \(f>l) <g> \<f) l p ) 



(20) 



for all i and icti- Each {T|^) (g> \<j)a)} contains No c- 
electrons and d-fermions with spin down and Nt — Nq 
c-electrons and d-fermions with spin up. The total num- 
ber of c-electrons and d-fermions in {T|^J,) ® \<t>a)} is 
Nt- Consider now the constraint of single occupancy of 
d-fermions. The number of c-electrons in this basis is 
N e = A A . The z-component of the total spin is 

S? ot T|^) ® = h Nt _ 2A o )T|0l) ® |^>. (21) 



Any lowest energy state of H' with 5f ot = \{N t — 2A ) 
can be expanded as 



i*(w)>=5Z w ^ T i^>®W>- 



(22) 



The coefficients {W 7 ^} can be regarded as a square 
matrix and chosen as hermitian since the transformed 
Hamiltonian TiJ'T^ possesses spin up-down symmetry. 

The proof of the positive definiteness of W is straight- 
forward, but a little tedious accordme. to the reflection 
positivity approach in the spin spacdj'Ej. The variational 
energy of H' is 

E{W) = {V(W)\H'\V(W)) 

= 2Tr(W 2 P) - U l Tr{WV l WV i ) 

(7 i£A 



£ \Ji\Tr(WV cdiai WV^ iai ) 

d,ai 

]T \Kij\Tr{WV diai j aj WVl ioiija .) 



i£f\d,oti 



MTr{WV diaa ,WV d 



^ XiSi(Si + 1), 



(23) 



where 



(P) a p = (CI E % c *> c .j> 

ijeA 



i£A d ,an 



ij£A d ,oti,oij 

E "^OWa - ^)Ka' i)CT - (24) 



and 
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(Vcdia^af) 
{V d lO-ijaj Jap 
{y^diaiQ.'. )a/3 



{r a \c ia di\^}, 

\d iryi _„d]^ T |<%). 



q | ^ICti . (7 



(25) 

(26) 
(27) 
(28) 



As W is hermitian and can be decomposed as W = 
V'DV where V^V = 1 and D is a diagonal matrix with 
the diagonal elements d a , denoted by D = diag{d a }. De- 
note \W\ = V^\D\V where \D\ = diag{\d a \}. \W\ is at 
least positive semidefinite as all its eigenvalues are non- 
negative. When all d a ^ 0, \W\ is positive definite. As 



Tr{W 2 P) = Tr(\W\ 2 P); 
Tr{WXWX^) < Tr(\W\X\W\X^) 

and all U > we have 

E{W) > E{\W\). 



(29) 
(30) 



(31) 



If |\t(W)) is the lowest energy state, l^dW])) must be 
also the lowest energy state in the variational principle. 
This indicates that one of the lowest energy states of H' if 
degenerate is at least positive semidefinite. Suppose \W\ 
is positive semi-definite (not positive definite) and there 
must exist one non-zero vector Vq such that |W|Vb = 0. 
From the Schrodinger equation for \W\ and using the 
variational principle, we have 



\W\N cl V = 0, if Ui > 0, 
\W\H V = 0, 

\W\V cdiai V = 0, ifJ^O, 

\w\v d V o =0, if Kij ^0, 
\W\V diaiai V =0, 



where 



(N ct ) a/3 = {4>a\n ia \4>p), 



(32) 
(33) 
(34) 
(35) 
(36) 



(37) 
(38) 



Eq. (36) is obtained from the single occupancy condition 
of d-fermions. When all Ui > and Ji ^ 0, the gen- 
eralised lattice is connected through t^ and Ji. In this 

we 
As 



case, we can show that for any non-zero vector V„,„ , v 
\W\V any 



by successively using Eqs. (|32 



have 

\W\ is an non-zero matrix, it is impossible that all of its 
eigenvalues are equal to zero. This indicates that \W\ is 
positive definite and non of its eigenvalues are equal to 
zero. If there are two positive definite states, for instance 
Wi and Wi , constructing the state W\ + cWi (c is a non- 
zero constant) and repeating the same procedure above 
we can show that W\ = W2 or W\ = — W2. Thus we con- 
clude that the lowest energy state in each subspace S^ ot 
is non-degenerate and positive definite. This conclusion 
holds even if some of J; = but the generalised lattice is 



connected through tij, Ji and Kij. However we cannot 
show the non-degeneracy of the lowest energy state in 
each subspace if the generalised lattice is not connected 
or not all U > except for some special cases. A de- 
tailed proof for the case of s = 1/2 is seen in Rcf. ||. In 
principle we have reduced the Kondo model with large 
spin ai to the model with spin 1/2. The proof for the 
case of 1/2 can be applied to the case of large spin. 

Since the lowest energy state is positive definite on 
the basis we choose., a theorem on a positive definite or 
semidefinite stateo'Ll can be applied to show Theorem (ii- 
iv). The theorem states: 

Given a positive semidefinite state |$) on a bipartite 
lattice with Na and Nb, then 

1) . If the state is an eigenstate of the total spin, the 
eigenvalue S is \N A - N B \/2 if \S Z \ < \N A - N B \/2, and 
\S Z \ otherwise; 

2) . The transverse spin-spin correlation function (for 
fermions with spin 1/2) obeys 

($|S+.S7|*) = e(<)e(j)Cy (39) 

where Cy > ( the equality only holds possibly in a 
positive semidefinite, not definite state). 

In our case, the difference of the total numbers of the 
two generalised sublattice A and B sites is J2ie/\ e (*) — 
Si6A d We get Theorem (ii) and (iv) combin- 

ing the theorem above and the positive definiteness of 
the lowest energy state. Theorem (iii) is also obtained as 
we have decomposed the spin Si into 2si 1/2-spins and 
the whole system is still on a generalised bipartite lattice. 
For example, 

<*|s+ • sji*) = J2 (*i4 i)T <w4a^.Ti*> 



<i)e{j)Fij (40) 



where 



Fij = Tr(W^V diaija .W(V diaija ^). (41) 



Fij > when the state is positive definite, and > when 
the state is positive semidefinite. 

In summary, we have provided several theorems on the 
ground state properties of the Kondo model at half-filling 
and in the case with large spin. The uniqueness and the 
total spin in the ground state are found. Furthermore 
we have also investigated the spin-spin correlation in the 
system. The co-existence of both antiferromagnetic and 
ferromagnetic long-range correlations is proved, which is 
similar to the case of spin 1/2. Hopefully, these exact re- 
sults are useful for understanding this highly correlated 
electron-impurity lattice. 
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